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@ An adaptive noise control system comprises a 
reference microphone (12) (Fig. 2) for generat- 
ing a reference signal (x(t)) that is con'elated 
with noise emanating from a primary noise 
source (10) , secondary loud spealcer sources 
(Si, 82, ... Sn) for generating a plurality of secon- 
dary sound waves, microphones (e^. 62 ... e^) for 
detecting a plurality of far-field sound waves in 
a farfield of the primary noise source and 
generating a plurality of error signals (e^ (t) , 
e2(t), ... eM (t)) each of which is indicative of the 
power of a corresponding far-field sound wave, 
and an adaptive controller (14) for controlling 
the secondary sources in accordance with the 
reference signal and the en-or signals so as to 
minimize the power in the far-field sound 
waves. 
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The present invention generally relates to the field of noise control and more particularly relates to adap- 
tive, active noise control systems. One preferred application of the invention is to control noise in a power gen- 
eration plant. ^ 
Free-field noise sources, such as internal combustion engines and combustion turbines, generate powerful 
5 low-frequency noise in the 31 Hz and 63 Hz octave bands (where the 31 Hz octave band extends from 22 Hz 
to 44 Hz and the 63 Hz octave band extends from 44 Hz to 88 Hz). Passive noise control requires the use of 
large, expensive silencers to absorb and block the noise. The size and cost of such silencers makes passive 
control unacceptable for many applications. An alternative to passive control is a combination of passive control 
and active control. Passive control abates noise better as the frequency of the noise increases and active con- 
ic? trol works better as the frequency of the noise decreases. Therefore a combination of passive and active control 
may advantageously be employed in many applications. 

The active control of sound or vibration involves the introduction of a number of controlled "secondary- 
sources driven such that the field of acoustic waves generated by these sources destructively interferes with 
the field generated by the original "primary" source. The extent to which such destructive interference is pos- 
ts sible depends on the geometric arrangement of the primary and secondary sources and on the spectrum of 
the field produced by the primary source. Considerable cancellation of the primary field can be achieved if 
the primary and secondary sources are positioned within a half-wavelength of each other at the frequency of 
interest 

One form of primary field that is of particular practical importance is that produced by rotating or recipro- 
20 eating machines. The waveform of the primary field generated by these machines is nearly periodic and, since 
it is generally possible to directly observe the action of the n^chine producing the original disturbance, the 
fundamental frequency of the excitation is generally known. Each secondary source can therefore be driven 
at a harmonic of the fundamental frequency by a controller that adjusts the amplitude and phase of a reference 
signal and uses the resulting 'filtered" reference signal to drive the secondary source. In addition, it is often 
25 desirable to make this controller adaptive, since the frequency and/or spatial distribution of the primary field 
may change with time and the controller must track this change. 

To construct a practical adaptive controller, a measurable error parameter must be defined and the con- 
troller must be capable of minimizing this parameter. One error parameter that can be directly measured is the 
sum of the squares of the outputs of a number of sensors. The signal processing problem in a system employing 
30 such an error parameter is to design an adaptive algorithm to minimize the sum of the squares of a number 
of sensor outputs by adjusting the magnitude and phases of the sinusoidal inputs to a number of secondary 
sources. S. J. Elliot et al.. in "A Multiple Error LMS Algorithm and Its Application to Active Control of Sound and 
Vibration." IEEE Trans, on Acoustics. Speech and Signal Processing, Vol. ASSP-35, No. 10, Oct 1987, de- 
scribe a least-mean-squares (LMS) based active noise control system, however that system converges too 
35 slowly for many applications. 

The present invention is directed to systenns for controlling both random and periodic noise in a single or 
multiple mode acoustic environment (In a multiple mode acoustic environment the amplitude of the sound va- 
ries in a plane perpendicular to the direction in which the sound propagates.) There are known systen^s for 
controlling random noise propagating in a single mode through a duct, however these systems do not work 
40 With multiple mode propagation. See U.S. Patents Nos. 4,044,203, 4,637,048 and 4,665,5498 and M. A. Swin- 
banks. The Active Control of Low Frequency Sound in a Gas Turbine Compressor Installation." Inter-Noise 
1982, San Francisco, CAMay 17-19, 1982.pp. 423-427. 

Accordingly, a primary goal of the present invention is to provide noise control methods and apparatus 
that can rapidly adapt, or converge, to an optimum state wherein the total noise received by a number of de- 
45 tectors placed in prescribed locations is minimized. Adaptive noise control systems in accordance with the pres- 
ent invention comprise reference means for generating a reference signal that Is correlated with noise eman- 
ating from a primary noise source, secondary source means for generating a plurality of secondary sound 
waves, detection means for detecting a plurality of far-field sound waves in a farfield of the primary noise 
source and generating a plurality of error signals each of which is indicative of the power of a corresponding 
50 far-field sound wave, and adaptive control means for controlling the secondary source means in accordance 
with the reference signal and the error signals so as to minimize the power in the far-field sound waves. 

In preferred embodiments of the invention, the reference means comprises means for detecting acoustic 
noise in the near-field of the primary noise source, the secondary source means comprises a plurality of loud 
speakers, and the detection means comprises a plurality of microphones. 
55 The adaptive control means in preferred embodiments comprises: (i) correlation means for generating au- 

tocorrelation data on the basis of the reference signal and generating crosscorrelation data on the basis of 
the reference signal and the error signals, (ii) FFT means for generating auto-spectrum data and cross- 
spectrum data on the basis of the autocorrelation and crosscorrelatron data, (iii) FIR means for filtering the 
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reference signal in accordance with a plurality of weighting functions and for providing filtered versions of the 
reference signal to control the output of the secondary source means, each weighting function being associated 
with a corresponding one of the secondary sound waves to be generated by the secondary source means, 
5 and (iv) adapting means for processing the auto-spectrum and cross-spectrum data so as to derive the weight- 
ing functions and for providing the weighting functions to the FIR filter means. 

Systems in accordance with the present invention may also advantageously comprise random number 
means for generating substantially random numbers and means for switching the input of the FIR means to 
the random number means. This enables the performance of a system identifi- cation function (described be- 
ta low) in accordance with the invention. 

The adapting means may comprise means for performing an inverse Fast Fourier Transformation of the 
said weighting functions prior to providing them to the FIR filter means. 

The present Invention may advantageously be applied in a power generation system comprising a conrv 
bustion turbine coupled to an exhaust stack. In such an application, an adaptive, active control system for cen- 
ts trolling multi-mode acoustic noise generated by the combustion turbine and emanating from the exhaust stack 
comprises reference means for generating a reference signal that Is correlated with noise generated by the 
combustion turbine, secondary source means for generating a plurality of secondary sound waves, detection 
means for detecting a plurality of far-field sound waves in a far-field of the exhaust stack and generating a 
plurality of error signals each of which is indicative of the power of a corresponding far-field sound wave, and 
20 adaptive control means for controlling the secondary source means in accordance with the reference and error 
signals so as to minimize the power in the far-field sound waves. 

The present invention also encompasses methods comprising steps corresponding to the respective func- 
tions of the elements described above. 

Noise control methods in accordance with the present invention can theoretically (i.e., under the right con- 
25 ditions) converge in one iteration. Moreover, systems in accordance with the invention are capable of efficiently 
achieving a large reduction in multi-mode noise, even in non-static noise environments. Other features and 
advantages of the invention are described below. 

Figure 1 is a schematic representation of a noise control system in accordance with the present invention. 
Figure 2 depicts a noise control system in accordance with the present invention in the context of a power 
30 generation system. 

Figure 3 is a more detailed block diagram of the noise control system of Fig. 1, with emphasis on the adap- 
tive control block 14. 

The theory underlying the present invention will now be described with reference to Figure 1 , which depicts 
a primary noise source NS surrounded by N secondary noise sources (or control sources) S^-Brj, where N rep- 

35 resents an integer. The primary noise source NS may be composed of one or more sources that radiate sound 
waves. Error microphones et-eM> where M represents a number greater than or equal to the number of sec- 
ondary sources N. detect sound waves in the far-field (approximately 1 50 ft. (45 meters)) of the prinnary noise 
source NS and provide feedback to a control system (not shown) that controls the secondary noise sources 
Si-Sn such that the total noise received by the error microphones is reduced. The secondary sources are driv- 

40 en by the output of a filter (not shown), which is part of the control system. The input to the filter, called the 
reference, may be derived by sampling the sound In the near-field of the primary noise source NS (e.g., within 
a few feet of NS). Alternatively, if the primary noise is periodic, a synchronization signal of a prescribed fre- 
quency may be used to generate the reference. Using current technology, the control system's filter can nrwst 
easily be implemented with a digital signal processor. The following analysis is therefore in the discrete time 

45 and "z" domains. (Those skilled in this art will recognize that the z domain is reached by performing a z-trans- 
form of sampled, or discrete time, data. The z-transformation of sampled data between the discrete time and 
z domains is analogous to the Laplace transformation of mathematical functions between the time and fre- 
quency domains. The z-transform Is a superclass of the discrete Fourier transform.) 

Referring to Figure 1 , an error microphone e^ (where m represents any number between 1 and M) receives 

50 sound from the primary noise source NS and the secondary sources S-t to Sn. The sound generated by NS 
and detected by error microphone e^ is represented as d^ in this analysis. Thus em(z) is given by the following 
equation: 
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Since there are M error microphones, the following matrix equation is formed: 

= + [y]S] (2) 

where, 

5 E] = [ei(z), e2(z). ... e^(z), ... eM(z)F 

D] = [d,(z). d2(z). ... d„(z). ...dM(z)r 
m = [YUz)] 

S) = [Si(z), S2(z), ... s„(z), ... Sn(z) r 
The element Y^n of the Y matrix represents the transfer function where the control signal Sn(z) is the input to 
10 secondary source Sn and the signal e^iz) is the output of error microphone e„,; i.e., Y^n = em(2)/S„(z) with 
Si(z). ... S„. i(z), S^i(z), ... Sn(z) = 0. 

As mentioned above, the control signal Sn(z) is the input to secondary source Sn, however it is also the 
output of the control system's digital filter (described below) with the input to the filter being a reference signal 
X(z). Sn(z) may be determined from X(z), a filter function Wn(z) and the following equation: 
^5 S„(2) = W„{z)X(z) (3) 

Substituting equation (3) into equation (2) yields: 

E] = D] + [Y]W]X(z) (4) 

where, 

W] = IWi(2). W2(z), ... Wn(z) F (5) 
20 The least squares solution to equation (4) (i.e., the values of W] that minimize the total noise power in E], 

given by ei2(2)+e22(z)+...eM2(z)) is 

W] = - ([Y]H[Y])-i[Y]HX • (z)DyX • (z)X(z) (6) 

where, 

[Y]" represents the conjugate transpose, or Hermitian, of [Y], and 
25 X*(z) represents the conjugate of X(z). 

In equation (6). the product X* (z)D] Is the cross-spectrum of the reference X(z) and the noise matrbc D]. 
The auto-spectrum X*(z)X(z) is a complex number and is divided into the cross-spectrum X*(z)D]. (Note that 
the cross- and auto-spectrums are also referred to in this specification as "Gxx(z)" and "GxemW. respectively.) 
The least-squares solution of W] can be found in one iteration with equation (6), provided there are no 
30 measurement errors in [Y], D] or X(z). In practice, however, errors in [Y], D] and X(z) are significant enough 
to require the following iterative solution: 

W] = W] . n([Y]H[Y^WX*(z)EJ/X•(z)X(z) (7) 
where ^ is a convergence factor. If ^ = 1 , equation (7) will reduce to equation (6) because E] = D] when W] = 
0]. Typical values of ^ are in the range of 0.1 to 0.5. 
35 Both the cross-spectrum X*E] and auto-spectrum X*X can be computed by taking the discrete Fourier 

transform, implemented, e.g., by the Fast Fourier Transform (FFT), of the crosscorretation of x(t) and e^W and 
autocorrelation of x(t), respectively (where x(t) represents the time-domain version of X(z)), The autocorrelation 
of x(t), designated Rxx(t), and crosscorrelation of x(t) and em(t). designated Rxem(t), are given by the following 
equations: 

40 



45 



L-l 



(8) 



50 



L-1 

^.e.(t) = lj^^e^{k)x{k-t) (9) 



where, 

55 k is the discrete time index, 

x(k) represents the reference signal in the discrete time-domain, 

em(k) represents the error signal, in the discrete time-domain, from error microphone number m, and 
L represents the number of samples used to compute Rxx(t) and Rxem(t) (note that the accuracy of the 
computation may be increased by increasing the number of samples L, however the disadvantage of making 
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L unnecessarily large is that the frequency at which the filters can be updated is Inversely proportional to L). 

To properly transform Rxx(t) and Rx©m(t) into the frequency donnain (i.e., the z-domain), the H-point vectors 
must be padded with zeros such that the resulting vector is 2H points long: 

5 

H zeros 



R«(t:) = [R«c(0), Rxx(l)/ Rxx(h), ... Rxx(H-l), 0, 0] 

Rxx(t) is then transformed to the auto-spectrum G^xCz) with a 2H-point FFT. Rxom(t) is transformed in the same 
manner to Gxem(z)- 

Due to causality constraints, the Wn(z) weighting functions must be transformed to the time-domain. The 
15 control signal Sn(t) is computed from 



- H-1 

S^{k) = v^(t)x{ic-t) (10) 

20 t=0 



where Wn(t) represents the time-domain versions of the filter functions Wn(z) and H represents the length of 
the fitter functions Wn(t) (also referred to as the number of taps in the respective filters). A 2H-point inverse 
25 discrete Fourier transform may be used to convert Wn(z) to Wn{t). Only the first H points of the result are used 
in equation (10). 

An application of the present invention to the suppression of noise emanating from the exhaust stack of 
a combustion turbine will now be described with reference to Figures 2 and 3. The dimensions of the cross- 
section of the stack are assumed to be greater than the wavelengths of the sound waves that emanate there- 
to from, therefore multi-mode noise will be generated. 

Figure 2 depicts a power generation system employing an active, adaptive noise control system in accor- 
dance with the present invention. In this system, a plurality of loudspeakers S^-Sn are positioned around the 
top rim of an exhaust stack 10 of a combustion turbine 11. A reference signal x(t) is measured by a probe mi- 
crophone 12 in the stack 10. A plurality of error microphones Bi-eM (with M >= N) are located in the far-field 
35 of the exhaust stack. An adaptive control system 14 takes feedback from the error microphones e^-eu and the 
reference signal x(t) from the probe microphone 12 and drives the loudspeakers S^-Sn so as to substantially 
cancel the noise detected by the error microphones. 

Figure 3 is a more detailed block diagram of the system of Figure 2, with emphasis given to the adaptive 
control system 14. (The turbine 11 and exhaust stack 10 are not shown in Fig. 3.) The reference numerals 12- 
40 42 refer to both structural elements (or hardware) and functional elements that may be implemented with hard- 
ware in combination with software; although the respective functional elements are depicted as separate 
blocks, it is understood that in practice more than one function may be performed by a given hardware element. 

The reference numerals are used as follows: 12-probe microphone, 14-adaptive control system. 16-switch, 
18-bus, 20-bus, 22-random number generator, 24-finite impulse response filters FIR^-FIRn, 26-secondary 
45 source loud speakers Si_Sn, 28-auto/cross-correlation blocks, 30-error detector microphones. 32-zero-pad 
blocks. 34-Fast Fourier Transform (FFT) blocks, 36-cross-spectrum array, 38-processing block, 40-processing 
block, and 42-inverse Fast Fourier Transform (IFFT) block. In one embodiment of the present invention, there 
are three processors (two digital signal processors and one microprocessor) involved in (1) filtering the refer- 
ence and generating the secondary source signals Si(t)-SN(t) (which drive the respective loudspeakers S^- 
50 Sn). (2) receiving the error signals and computing the autocorrelation and crosscorrelation vectors Rxx(t), 
F^xei(t)-RxeM(t), and (3) carrying out the FFTs, updating the filter coefficients and carrying out the inverse FFT. 

One problem encountered by the present inventors is the causality of the reference signal with respect to 
the sound at the secondary sources. The group delay characteristics of the low-pass fitters (LPFs), the high- 
pass response of the secondary sources S^-Sn. and the delay of the digital filters FIRi-FIRn must be less than 
55 the time that the noise takes to travel from the probe microphone 12 to the closest secondary source. There- 
fore, to derive each control signal Sn(t) the reference signal x(t) is filtered in the time domain with a finite impulse 
response (FIR) filter. It has been argued that the filters are best implemented by an infinite impulse response 
(MR) filter. See L. J. Eriksson, et al., The Selection and Application of an MR Adaptive Filter for Use In Active 
Sound Attenuation," IEEE Trans, on Acoustics, Speech and Signal Processing, Vol. ASSP-35, No. 4, April. 
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1987, pp. 433-437 and L J. Eriksson, et al-. The Use of Active Noise Control for Industrial Fan Noise," American 
Society of Mechanical Engineers Winter Annual Meeting, Nov. 27 - Dec. 2, 1988, 88-WA/NCA-4. However, be- 
cause of the potential instability of IIR filters^ the present invention employs intrinsically stable FIR filters, with 

5 the understanding that a large number of filter taps may be required in particular applications. 

Another problem encountered by the inventors is the updating of the filter coefficients Wr,(t) of the FIR 
filters. Typically, adaptive filters implemented in the time-domain are updated in accordance with time-domain 
algorithms. Elliot describes such a system in S. J. Elliot, et al., "A Multiple Error LMS Algorithm and Its Appli- 
cation to Active Control of Sound and Vibration," IEEE Trans, on Acoustics, Speech and Signal Processing, 

10 Vol. ASSP-35, No. 10, Oct. 1987. However, the convergence time of an LMS-based control system the time 
that the control system 14 needs to adjust the filter coefficients to optimum values) can be many orders of 
magnitude greater than the convergence time of the present invention, which adjusts the filter coefficients in 
the frequency domain. 

Frequency domain adaptive algorithms have very advantageous properties, such as orthogonal reference 
15 signal values, which are a direct result of taking the FFT of the autocorrelation of x(t) (i.e., the frequency com- 
ponents of Gxx(z) are independent of one another). In addition, the entire updating process is decomposed into 
harmonics, or frequency "bins", which makes the process easier to understand, and thus control, than a time- 
domain process. In preferred embodiments of the present invention, the filter functions Wi(z)-W|si(z) are gen- 
erated in the frequency domain and then converted to the time-domain functions Wi(t)-WN(t). The time-domain 
20 functions Wi(t>-WN(t) are provided via a set of busses 20 (only one bus 20 is shown in Fig. 3) to the FIR filters 
FIRi-FIRn. 

The adaptive control system 14 mustf irst /den^^y the system before optimizing the FIR filters. System iden- 
tification involves determining the respective transfer functions Ymn(t) from the inputs of the digital-to-analog 
converters (DACs) (Fig. 3), through the speakers Sh, the acoustic path from the speakers Sn to the error mi- 

25 crophone 6^, and finally to the outputs of the analog-to-dtgital convertor (ADCs). This is accomplished by gen- 
erating random numbers with a digital random number generator 22 and outputUng these numbers via a switch 
18 to a bus 18 coupled to the respective FIR filters and to inputs of autocorrelation and crosscorrelatkin blocks, 
which compute autocorrelation and crosscorrelation data. As a final step, the auto- and crosscorrelation data 
(Rxx(t) and Rxei(t) -RoxmCO) is converted to 2H-point frequency-domain data (Gxx(z) and Gxei(z)-Gx©M(z)) by zero- 

30 pad and FFT blocks 32, 34. 



System Identification 



The system identification process may be summarized as follows: 
35 Step 1 : Set switch 16 to the random number generator 22. 
Step 2: Setn=1 

Step 3: Zero all FIR coefficients W„(h) (for h = 1 to 2H-1) and set Wn(0) to 1 .0. 
Step 4: Compute autocorrelation and crosscorrelation data using equations (8) and (9). 
Step 5: Zero pad Rxx(t) and Rxoi(t)-RxeM{t) and take the FFT of each to produce Gxx(h) and Gxei(h)- 
^ GxoM(h), where h now represents the harmonic index of the FFT and takes values from O to 2H- 

1 . (Note that the actual frequency corresponding to the index h is a function of the sampling fre- 
quency and the number of points 2H, and may be determined t>y well-known techniques.) 
Step 6: Compute Ymn (h) using the following formula: 

YUh) = GxeJh)/Gxx(h) (11) 
45 for h = 0 to H-1 and m = 1 to M. 

Step 7: If n is not equal to N (the number of secondary sources), increment n by 1 and repeat steps 3 
through 6. 

Step 8: Compute the Z matrix for each harmonic h as follows: If N = M. compute Z as follows: 
For h = 0 to H-1 do 
50 IZ^„(h)l = [Y^„(h)hi (12) 



If M > N, compute Z as follows: 
For h = 0 to H-1 do 

[Z„,„(h)l = [Y„„(h)]HIY(h)hi[Y(h)]H (13) 
55 (Note that the superscript H In equation (13) represents the Hermitian operator.) 



Adaptation 

Adaptation determines the optimum filter coefficients for each FIR filter. The adaptation process may be 
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15 



summarized as follows: 
Step 1 
Step 2: 
Step 3 
Step 4; 



Step 5: 



Step 6: 
Step 7: 
Step 8: 



Set switch 16 (Fig. 3) to the ADC of the reference channel coupled to the probe microphone 12. 
Zero all FIR coefficients Wn(t) and Wn(z) forn = 1 to N. 

Compute autocorrelation and crosscorrelation data using equations (8) and (9). 
Zero pad Rxx{t) and Rxei(t)-RxeM(t) to 2H points and take the FFT of each to produce Gxx(h) and 
Gxei(h)-GxeM(h); setn = 1. 

Compute frequency-domain filter coefficients Wn(h)] using 

Wn(h)l = W„(h)] - ^[Z(h)]G^e(h)]/(G«Mh)G«(h)) (14) 
where h =0to2H-1. 

Inverse discrete Fourier transform Wn(h)] into the time-domain coefficients W„(t)]. 
Load updated time-domain coefficients Wn(t) into filter FIRn. 

If n is not equal to N (the number of secondary sources), increment n by 1 and repeat steps 5 
through 7. 



Necessary Conditions for Active Control 



The following conditions must be met for active control to successfully reduce random noise (these are 
designated the %ur C's"): 

20 1) There must be sufficient coherence between the reference microphone signal and the far-field sound 

pressure. 

2) If there are multiple noise sources, they must have coalesced and appear as one source. 
3} Sampling of the reference signal must be sufficiently advanced in time to compensate for the transient 
response of the active control system. This is called the causality requirement. 
25 4) The secondary control sources must have sufficient capacity to generate a cancelling sound field. 

Each of these requirements are briefly discussed below. 



Coherence 



30 The coherence between two signals ranges from 0 to 100 percent. In the case of the exhaust stack 10, 

the reference microphone 12 detects the sound inside the stack and, barring any other noise sources, this 
sound should be highly related to, or coherent with, the sound at the top of the stack and the sound detected 
by the far-field microphones Si-Om. In other words, the sound power detected by the far-field microphones 
should nearly be 100% the result of the sound radiating from the top of the exhaust stack 10. In reality, however, 

35 the percentage of the sound power detected by the far-field microphones that comes from the top of the stack 
drops as the sound generated by other unrelated noise sources (such as a mechanical package, turbine inlet 
and turbine housing) is detected. For example, if the coherence between the sound at the top of the stack 10 
and the sound detected by the far-field microphones Oi-Om is 60%, then 40% of the sound power detected in 
the far-field will be related to other noise sources, such as the turbine housing and mechanical package. To 

40 illustrate the importance of coherence in assessing the value of a given noise control system, suppose that 
all of the noise radiating from the exhaust stack were eliminated. Then the sound power in the far-field would 
decrease by 60%, or 4 dB. 

The following table lists the theoretical maximum noise reduction for a given coherence between the ref- 
erence signal x(t) and the far-field signals. 



45 



50 



55 



COHERENCE 


NOISE REDUCTION POWER RATIO 


100% 


Infinite 


99% 


20 dB 


90% 


10 dB 


80% 


7dB 


60% 


4dB 


50% 


3dB 
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Coalescence 

The combustion chambers of a combustion turbine can be considered distinct and mutualiy incoherent 
5 noise sources. The sound emanating from each of the combustion chambers mixes, or coalesces, as it prop- 
agates through the exhaust section and into the exhaust stack. Once the noise has coalesced in the exhaust 
stack, the sound at any location in the stack should be more than 90% coherent with the sound at any other 
location in the stack. However, turbulence noise produced, e.g., by the flow of exhaust gases through the ple- 
num and silencer creates spatially incoherent noise In the exhaust stack and thus the coherence between the 
10 sound at two points in the stack will decrease as the distance between the two points increases. Turbulence 
noise generated by flow through a silencer is often called self noise. If the exhaust flow is turbulence-free after 
the exhaust silencer, the spatially incoherent sound at the exhaust silencer will coalesce once again as it prop- 
agates up the exhaust stack. 

IS Causality 

Causality refers to the requirement that the reference signal x(t) must be obtained a sufficient amount of 
time before the sound reaches the control speakers S^Su for the control system 14 to filter the reference signal 
and drive the speakers. The transient delay of one embodiment of the control system is about 3 milliseconds 
20 (ms) and the speakers have a transient delay of at>out 12 ms; Therefore the total time delay from the reference 
microphone input to the acoustic output of the speakers is about 15 ms. Since sound travels about 1 foot per 
1 ms, the reference microphone should be approximately 15 ft (4.6 meters) from the top of the stack. Ashorter 
distance may produce satisfactory results for some applications. 

25 Capacity (or Control Power) 

The loudspeakers S^-Sn, should be able to generate as much sound power as that emanating frmi the 
stack 10. However, because of the interaction between independent control sources, the specified power levels 
for the loud-speakers should be at least twice that radiated by the exhaust stack. 

30 

Experimental Results 

A three speaker (i.e., three secondary sources) and four error microphone active control system in accor- 
dance with the present invention has been tested. Alow-pass-filtered (0-100 Hz) random signal acted as the 
35 driving signal to a primary noise source speaker and as the reference signal x(t). The filter coefficient optinv 
ization process was frequency-limited by the operator to 20-170 Hz. Reductions in sound pressure level (SPL) 
of up to 27 dB were achieved between 20 Hz to 120 Hz. A slight increase in SPL was noted between 120 Hz 
and 160 H2l This problem was solved by setting the upper frequency limit to 120 Hz. 

40 

Claims 

1. A power generation system, comprising a combustion turbine (11) coupled to an exhaust stack (10), and 
an adaptive, active control system for controlling multi-mode acoustic noise generated by said combustion 
45 turbine and emanating from said exhaust stack, said active control system characterized by: 

(a) reference means (12) for generating a reference signal (x(t)) that is correlated with noise generated 
by said combustion turbine; 

(b) secondary source means (Si, Sa, ... S^) for generating a plurality of secondary sound waves; 

(c) detection means (Si, 02, ... Om) for detecting a plurality of far-field sound waves in a far-field of said 
50 exhaust stack, and generating a plurality of error signals (ei(t), e2(t), ... eM(t)) each of which is indicative 

of the power of a corresponding far-field sound wave; and 

(d) adaptive control means for controlling said secondary source means in accordance with said ref- 
erence signal and sakJ error signals so as to minimize the power in said far-field sound waves, saM 
adaptive control means comprising: 

55 (i) correlation means for generating autocorrelation data (Rxx(t)) on the basis of said reference signal 

and generating crosscorrelation data (Rxei(t). Rxe2(t). RxeM(t)) on the basis of said reference signal and 
saki error signals; 

(ii) FFT means for generating auto-spectrum data (Gxx(h)) and cross-spectrum data Gxei(h), Gxe2(h), ... 
GxeM(h)) on the basis of said autocorrelation and crosscorrelation data; 
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(iti) FIR means, coupled to said reference means, for filtering said reference signal in accordance with 
a plurality of weighting functions (wi(t), W2(t). ... WnW) and for providing filtered versions of said refer- 
ence signal to control the output of said secondary source means, each weighting function being as- 
sociated with a corresponding one of said secondary sound waves to be generated by said secondary 
source means; and 

(iv) adapting means for processing said auto-spectrum and cross-spectrum data so as to derive said 
weighting functions, and for providing said weighting functions to said FIR means, 
(d) adaptive control means for controlling said secondary source means in accordance with said ref- 
erence signal and said error signals so as to minimize the power in said far-field sound waves. 

A power generation system as described in claim 1, further characterized in that said reference means 
comprises means for detecting acoustic noise in the nearfield of said exhaust stack. 

A power generation system as described in daim 2, further characterized in that said secondary source 
means comprises a plurality of loudspeakers. 

. A power generation system as described in claim 3, further characterized in that said detection means 
comprises a plurality of microphones disposed In the farfield of said exhaust stack. 

A power generation system as described in dalm 4, further characterized in that random number means 
is provided for generating substantially random numbers and means for switching the input of said FIR 
means to said random number means, whereby a system identification function may be performed. 

L A power generation system as described in daim 5, further characterized in that said adapting means 
further comprises inverse FFT means for performing an inverse- Fast Fourier Transformation of said 
weighting functions priorto providing them to said FIR means^ 

7. A method for controlling noise emanating from a primary noise source, said method characterized by 
the steps of: 

(a) generating a reference signal that is correlated with noise emanating from said primary noise 
source; 

(b) generating a plurality of secondary sound waves in a near-field of said primary noise source; 

(c) detecting a plurality of far-field sound waves in a far-field of said primary noise source, and gener- 
ating a plurality of error signals each of which is indicative of the power of a corresponding far-field 
sound wave; and 

(d) controlling the generation of said secondary sound waves in accordance with said reference signal 
and said error signals so as to minimize the power in said farfield sound waves, said controlling step 
induding the following sub-steps: 

(1) generating autocorrelation data on the basis of said reference signal and generating crosscor relation 
data on the basis of said reference signal and said error signals; 

(ii) generating auto-spectrum data and cross-spectrum data on the basis of said autocorrelation and 
crosscorrelation data; 

(ill) processing said auto-spectrum and cross-spectrum data so as to derive a plurality of weighting 
functions; and 

(iv) filtering said reference signal in accordance with said weighting functions, and employing filtered ver- 
sions of said reference signal to control the generation of said secondary sound waves, each weighting func- 
tion being associated with a corresponding one of said secondary sound waves to be generated. 

8- A method as described in daim 7, further characterized in that step (a) comprises detecting acoustic noise 
in the near-field of said primary noise source. 

9. A method as described in claim 8. further characterized in that step (b) comprises the excitation of a plur- 
ality of loudspeakers. 

10. A method as described in daim 9, further characterized in that step (c) comprises the detection of said 
far-field sound waves with a plurality of microphones disposed in the far-field of said primary noise source. 

11. A method as described in daim 10, further characterized in that said adapting step (d) (iv) comprises per- 
forming an Inverse fast Fourier transformation of said weighting functions. 
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